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ABSTRACT 

We present a kinematic analysis of the dense molecular gas in the central 200 parsecs of the nearby 
galaxy NGC 1097, based on Cycle observations with the Atacama Large Millimeter/sub-millimeter 
Array (ALMA). We use the HCN(4-3) line to trace the densest interstellar molecular gas («h 2 ~ 10 8 
cm~ 3 ), and quantify its kinematics, and estimate an inflow rate for the molecular gas. We find a 
striking similarity between the ALMA kinematic data and the analytic spiral inflow model that we 
have previously constructed based on ionized gas velocity fields on larger scales. We are able to follow 
dense gas streaming down to 40 pc distance from the supermassive black hole in this Seyfert 1 galaxy. 
In order to fulfill marginal stability, we deduce that the dense gas is confined to a very thin disc, and 
we derive a dense gas inflow rate of 0.09 M Q yr~ 1 at 40 pc radius. Combined with previous values 
from the Ha and CO gas, we calculate a combined molecular and ionized gas inflow rate of ~ 0.2 
M yr _1 at 40 pc distance from the central supermassive black hole of NGC 1097. 
Subject headings: galaxies: kinematics and dynamics — galaxies: active — galaxies: individual 
(NGC 1097) 



1. INTRODUCTION 

The central region of the nearby barred Seyfert 1 
galaxy NGC 1097 displays a number of intriguing mor- 
phological and kinematic features. At ~ 1 kpc, an al- 
most circular ring-like feature marks the transition be- 
tween the prominent R ~ 8 kpc galactic bar and the rel- 
atively diffuse region interior to the ring. The bar hosts 
two prominent dust lanes, bot h originating at around th e 



dKohnoet all2003; Hsieh et aO 2ljM[20lT] [20121. More 



Corotation radius of the bar (Pinol-Ferrer et al. 2013) 



cutting through the inner ring and transforming into nu- 
clear spirals that co ntinue down to ~ 3.5 pc distance from 



over, Pinol-Fe rrer et al.| ( |201ip showed that the bulk of 
the interstellar gas at the centre of NGC 1097 (like in 
many other galaxies) is in the molecular phase, and there- 
fore, a detailed analysis of the central gas concentration 
needs to account for the different physical conditions. 
The interplay between the different phases provide cru- 
cial clues to understanding the energies involved in re- 
distributing the gas in a way that leads to the observed 
phase transition efficiencies. 
The discovery of bro ad (~ 10000 km/s) double-peaked 



the active nucleus (Lou et al. 2001 Fathi et al. 2006) 



Ha emission lines by Storchi-Bergmann et al. 



The dust lanes are accompanied by diffuse ionized ga s ™ ak f NGC 1097 also an ideal laboratory tor stu^yVg 



(1993) 



revealing clear kinematic signatures of bar-induced gas 
inf low ov er the entire face of the galaxy ( Pinol-Ferrer et 
aL] [2013l ). 
N eutral gas and i onized gas data cubes ( Ondrechen et 



al.|1989 



Fathi et al.|2006[ ) have confirmed the 'abundant 
presence of these two phases of the interstellar medium 
across the central kpc radius. However, interferomet- 
ric molecular gas maps show emission confined to the 
nuclear ring and the central 2-3 hundred parsecs radius 



the fate of the gas accumulated in the centers of active 
galactic nuclei (AGN). At a distance of 14.5 Mpc (i.e., 
~ 70 pc/"), this galaxy is also suitable for high-resolution 
studies of the physical processes that cause the mate- 
rial/fuel to loo se its angular momentum and fall toward 
the AGN (e.g., |Storchi-Bergmann et al.|2003 |. 

Although it is straightforward to transport gas down 
to the central kpc and induce enhanced star formation, 
it is more difficult to make the gas reach smaller scales 



Fat hi et al. 



(few pc) required to fuel an AGN. In rotating systems, 
perturbations can cause the potential to become non- 
axisymmetric, and torques exerted by the subsequent 
non-axisymmetr ic features are a ble to drive material to- 
ward the centre ( |Schwaxz|1984| . |Sheth fc Teubenj |1999| ) 
have argued that magnetic stress may aid the infalling 
gas to complete the last few parsecs down to the central 
supermassive black hole (SMBH). 

To build a realistic scenario for the fate of the gas that 
is piling up around an AGN to eventually fuel it, one 
has to make a detailed analysis of the distribution and 
kinematics of multiple phases of the interstellar gas in 
the region of interest. 

Based on ionized gas kinematic maps (van de Ven 
& Fathi 2010, hereafter vdVF, and Piiiol-Ferrer et al. 
2013), we have derived a concise picture for NGC 1097, 
according to which, the gravitational perturbation that 
once gave rise to the formation of the prominent bar 
drives the evolution of structure inside the bar as well 
as the outer spiral arms. The outer spiral arms are con- 
fined between the Corotation radius and reach beyond 
the Outer Lindblad Resonance radius of the main galac- 
tic bar. The circumnuclear ring once forme d at the loca- 
tion of the O uter Inner Lindblad resonance ( Pihol- Ferrer 
et al. 2013), and has likely migrated toward the centre 
of the^alactic gravitational poten tial (|Regan fc Teuben 



2003|lvan de Ven fc Chang||2009[ ). Inside the ring, the 
non-circular velocities are consistent with the presence 
of two spiral arms (in morphology). The nuclear arms 
are well disguised in optical images, and several image- 
enhancement techniques have led differe nt authors to ar- 



gue for a different n umber of arms (e.g., Lou et aL]|2001 



Davies et al. 2009 and vdVF). Most of the disagree 
ments concern the inner ~ 100 pc, where increased dust 
content may have distorted both images and kinematic 
measurements in the optical and near infrared. 

Here we present a quantitative analysis of the kine- 
matics of the densest molecular gas within the cen- 
tral kpc radius of NGC 1097 based on Atacama Large 
Millimeter/sub-millimeter Array observations of Hydro- 
gen cyanide, HCN(4-3) (Proposal 2011.0.00108. S, PI: 
Kohno). We compare the ALMA kinematic data with 
a dynamical model that we have previously constructed 
based on two-dimensional spectroscopic data of ionized 
gas. 

2. THE DATA 

NGC 1097 was observed with ALMA on November 5 
and 6, 2011, with 14 and 15 antennas, respectively. Our 
Band 7 observations targeted the HCN(4-3) line at rest 
frequency of 354.505 GHz with the original channel spac- 
ing of 488.28125 kHz. To increase the signal strength, the 
data were binned by a factor 20, leading to an effective 
channel width of 9.77 MHz (~ 8.3 km s _1 ). The primary 
beam was 18" 1 with a synthesized beam of l'/50 x 1'.'20 ( 
~ 105 x 84 pc) and at —72.4° position angle, sam pled at 
O^/p i x. Al l data specifications are described in Izumi 



et al.| ( |2013[ ). 

Figure ffihows the integrated intensity map of HCN(4- 
3) with the velocity moment maps and models, described 
below. At this resolution the spiral arms are difficult 
to see in the integrated intensity map, which sums all 
velocities, but are easier to detect in the kinematics. 

The HCN(4-3) kinematics were derived in two ways. 



Simple first Moment map (intensity-weighted mean ve- 
locity map) was cross checked with Gaussian fitting to 
each individual spectrum. We found no significant sig- 
natures of non-Gaussianity, however, the Gaussian fits 
resulted in a generally noisier velocity field. We use the 
Moment 1 maps and exclude all pixels for which the cor- 
responding spectrum Amplitude-over-Noise A/N < 20 
(see Fig. nj). We find almost no signal outside the area 
illustratea in Fig. fl] at lower A/N values. This confirms 
that the dense gas is co nfined to a small region around 
the AGN in NGC 1097 ( THsieh et al.|[2012] ). 

The ratio between the Einstein coefficient and colli- 
sional rates for the HCN(4-3) transition results in an 
estimated critical density of a few times 10 7 to 10 8 cm -3 
with a small dependency on temperature. Our calcula- 
tion is conformal with previously reported densit ies of 
nuo ~ 10 8 cm~ 3 at 40 K kinetic te mperature (e.g., Choi 

Such high densities 



et al. 



are a 



2000 Takakuwa et al. 2007) 



so consistent with the deep o bscuration expected 



toward the center of a galaxy (e.g., |Kohno et aL||1996 
Sakamoto et al.|2010[ ). 



lb constrain the large scale gas kinematics, we have 
used a mosaic of Fabry-Perot interferometric observa- 
tions at / /83 spatial sampling covering a 7' x 7' field at 
15 km s -1 spectra l resolution (Dicaire et al.|2008 Pifiol 



Ferrer et al.||2013[ ). To further look into the central few 
100 parsecs, we used Ha two-dimensional velocity field, 
at CK'l spatial sampling and 85 km s" 1 spectral resolu- 
tion, from the Gemini South Te lescope's Integral Field 
Unit, covering the inner 7" x 15" ( jFathi et al.|2006| ). The 
combination of the two sets of data is imperative, as they 
present a coherent dynamical model for NGC 1097 from 
~ 20 kpc down to ~ 100 pc from its central SMBH. 

3. KINEMATIC PARAMETERS 

We assume that the galactic disc is predominantly ro- 
tating and apply the prescription used in vdVF to quan- 
tify the velocity field shown in the top middle panel of 
Fig. [T] We divide the observed velocity field into con- 
centric rings, each containing > 7 pixels (this sets the 
limit for the innermost radius to 0'.'55 < 40 pc). We fix 
the inclination at 35° and apply a x 2 minimization to 
obtain the central coordinates (xo,J/o)j systemic veloc- 
ity V^y S and position angle PA of the disc. After fixing 
these parameters, we make the final \ 2 fit to the desired 
modes (here up to and including 3rd order) of the Fourier 
decomposition mathematically formulated as 

Vios = Kys + ^2,[c n {r) cosnd + s n (r) sinn6>] sini. (1) 
n=l 

We apply Eq. ^ to the HCN(4-3) velocity field and 
find an agreementbetween the kinematic and photomet- 
ric center to within one pixel (see Fig. IT]). Each ring is 
populated at more than 225 degrees, hence, we ensure 
that we do not need to assume any level of symmetry 
in the observed non-circular motions. We derive a sys- 
temic velocity of 1300 ± 48 km s _1 and the kinematic 
position angle of 147 ±6°. In Fig. [TJwe illustrate the dif- 
ferent stages of the fits described here, and a three-fold 
symmetry can be found in the non-circular velocity field 
(as predicted for this region in NGC 1097 by vdVF). The 
avreage final residual velocities are 10 km s" 1 , indicat- 
ing that we have reproduced 95% of the observed velocity 
features with the Fourier decomposition method. 
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Fig. 1. — Left: Hubble Space Telescope Ha image with a structure map showing the nuclear dust features | |Fathi et al.||2006y and the 
footprint of the ALMA HCN(4-3) pixels with A/N > 20 (magenta contours). Overplotted, in green, are the spiral arms modeled by 
vdVF. Right: ALMA HCN(4-3) distribution and velocity field. The circular and non-circular motion maps, Fourier reconstructed and final 
residual maps are displayed in km s _1 . The black curves in the bottom left panel outline the predicted location of the non-circular motions 
due to the model spiral arms corresponding to the two-arm morphological arms drawn on the HST image (in green). The non-circular 
velocities follow a similar three-fold symmetry. The ellipse at the top left panel displays the beam size, and the cross marks the kinematic 
centre, offset from the peak centroid by less than 0'.'3. 

Fitting Gaussians to each spectrum gives us the aver- 
age velocity uncertainty for all the pixels at 15 km s _1 . 
These errors are used to derive the uncertainties for the 
derived kinematics by means of Monte Carlo simulations. 
Repeated application to the Gaussian-randomized veloc- 
ity field yields the uncertainties on the Fourier param- 
eters. Our simulations show th at the average ro tation 
curve uncertainty is 20% (see also Fathi et al. 2005 ). Sim- 
ilarly, the higher Fourier term uncertainties have been 
calculated, and we plot the first, second and third Fourier 
terms of the HCN(4-3) data together with those derived 
from the GMOS data in Fig. [2j 



4. ANALYSIS 

The ionized gas velocity Fourier terms were modeled 
by vdVF, who constrained the detailed structure of the 
nuclear spiral arms and the associated gas inflow from 
kinematic data (solid curves in the lower panels of Fig.[2]). 
Their nuclear spiral structure is consistent with a weak 
perturbation in the gravitational potential due to a two- 
arm logarithmic spiral (in morphology) with a pitch angle 
of 52 ± 4° derived directly from the Fourier expansion of 
the model velocity field. Similarly large pitch angles in 
the very central parts of galaxies have also been modeled 



by |Yuan fc Yang (J2006J). Furthermore, the innermost 
~ 100 pc radius of the data points analyzed by vdVF 
displayed the largest errors in the third Fourier terms, as 
a consequence of dust contamination. 

We test the effect of beam smearing on the derived 
kinematic parameters by artificially smoothing by a fac- 
tor two, and find that the Fourier terms remain virtually 
unchanged. A notable effect of beam smearing is that 
it may lead to incorrect kinematic center, which in turn 



values would be affected by beam smearing. 

In Fig. [2j we overplot the ALMA HCN(4-3) velocity 
Fourier components on those derived from the GMOS 
data. The top panel reveals a very good agreement be- 
tween the rotation curve (~ c\ in Eq 1) of the HCN(4-3) 
and that of the Ha gas in the central 200 pc. Hence, 
due to the agreement of their overall kinematics, the 
bulk rotation of the dense gas is co-planar to the ionized 
interstellar gas. Further light can be shed on the dy- 
namical behaviour of the HCNY4-3) by considering the 
higher Fourier terms. Figure. [2] illustrates the analytic 
spiral model that has lead to the kinematic derivation 
of the pitch angle of the nu clear spirals with associa ted 
gas inflow rates (vdVF and Pinol- Ferrer et al. 2011). A 
simple overplot of the ALMA data (filled circles) dis- 
plays the striking agreement between the relatively non- 
contaminated HCN(4-3) kinematic behaviour and the 
model predictions. The analytic spiral model predicts the 
HCN(4~3) velocity Fourier terms down to CK'55 < 40 pc. 

Besides the n — 1 and n = 3 Fourier terms in the veloc- 
ity field, the non-circular motions also contain marginal 
n = 2 terms. vdVF argued that the n = 2 terms in the 
Ha were most likely due to dust contamination and pos- 
sible shocks associated with the gas streaming along the 
nuclear spiral arms (see also Valotto & Giovanelli 2004). 



may c ause uncertainties in the third Fourier terms. Wong 
(2000) found that smearing of 10"could produce up to 10 
km s _1 third Fourier terms. This is less than our error 
bars. In light o f these tests and following the discussion 



in (Wong 2000 Chapter 2), it unlikely that the innermost 



This seems to be confirmed here with the HCN(4-3) sec- 
ond terms consistent with zero inside 100 pc radius. 

5. DENSE GAS STREAMING DOWN TO < 40 PC FROM 
THE SMBH 

We use the derived C\ curve to measure the dynamical 
mass, assuming a thin disc model and velocity uncer- 
tainties of 20%. We derive the mass within the central 
100 pc to 3.5+2:o x 1q8 m ©' and inside the 40 P c radius, 
8. 0+3'g x 10 6 M . These masses are based on velocity 
measurements well outside the s phere of influence of the 
SMBH (~ 13 pc, |Peebles||1972[ ). Even so, in theory, a 
significant contribution from the SMBH is expected on 
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Fig. 2. — Comparing the quantified kinematics of the ionized 
gas in the optical (crosses) and the HCN(4-3) ALMA data (filled 
circles) to the analytic model prediction presented in vdVF. The 
spiral model is not a fit to any of the two data sets shown here, 
but a carefully chosen set of spiral parameters that best match 
the GMOS data. The HCN(4-3) data points are simple overplots 
of the new ALMA data showing a striking agreement with the 
model predictions. The model does not contain C2 or S2 terms (see 
Eqs (A12) and (A13) in vdVF). 



the rotating velocities at 40 pc radius. However, the 
synthesized beam size of the ALMA data almost entirely 
smears out the effect of the SMBH at this radius, and 
therefore, a rotation curve at higher spatial resolution 
is needed to measure the dynamical mass of the SMBH. 
Furthermore, it should be noted that our mass estimate 
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Fig. 3. — Disc scale height over radius h/R as a function of 
radius for the HCN(4-3) emitting gas at nj] 2 ~ 10* (solid curve). 
For comparison, we also illustrate this parameter for the ionized gas 
(dotted curve) and the molecular counterparts at lower densities 
(dashed curve). 



from the observed cold molecular gas rotation curve are 
not corrected for the contribution of non-cirular motions 
in the ci curve (see Fig. 3 in vdVF). Thus, these are 
likely lower mass limits. 

The prescription for deriving mass inflow rates along 
the nuclear spirals has been presented in vdVF. They 
combined the resulting inflow velocity corresponding to 
the model spiral arm parameters, with the gas density 
in the spiral arms inferred from [Sn] emission line ratios. 
They calculated the ionized gas inflow rate as a function 
of radius, reaching 0.033 M Q yr _1 at a distance of 100 pc 
from the central SMBH. The inflow rate was later refined 



by Pinol- Ferrer et al. ( 20TT| ) who used CO gas measure- 
ments to derive one order of magnitude higher molecular 
gas inflow rate for a margina lly stable disc model. Using 
this formalism, IPihol-Ferrer et all (120111) derived the CO 



gas inflow of 0.3M Q yr _i at 100 pc radius 

As we compare the HCN(4-3) kinematics with our pre- 
vious resul ts for ionized and CO gas (vdVF, Pinol-Ferrer 
et al.|2011 |, we set the sound speed at 10 km s _i and de- 
rive the dense gas scale height (Fig. [3]) using the observed 
epicyclic frequency. We then use the kinematic parame- 
ters from the analytic spiral model to trace the gas inflow 
down to the resolution limit (0 / .'55 < 40 pc). The data 
set at hand does not allow investigating density differ- 
ence between the nuclear spiral arms and the inter-arm 
region in NGC 1097. We assume 10% overdensity, simi- 
lar to the upper limit of the observed value from the [Sn] 
doublet, and in agreement w ith numerical simulations 
( |Englmaier fc Shlosman|2000| . We derive the HCN(4-3) 
inflow rate of 0.3 M yr _1 at 100 pc, and 0.09 M yr" 1 
at 40 pc from the SMBH ( see Fi g. |4| ■ The marginal sta- 
bility criterion of Romeo ( 1994 1 ensures that if the gas 
reaches higher densities, it will be confined to a thinner 
disc. Hence, changing the gas density by two orders of 
magnitude will not change the mass inflow rate (Fig.pl. 
However, accounting for the minimum density variation 
(arm versus inter-arm) derived from the [Sn] line ratios 
(i.e., 5p = 5%), the inflow rate would decreased by a 
factor 2 (Fig. [2]). 

The critical value 0.01 MEdd for the tr ansition be tween 
FINER and Syl galaxies was found by Ho (2005) with 
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Fig. 4. — Dense gas inflow rate derived from the analytic spiral 
model of vdVF as a function of galactocentric radius. Changing 
the molecular gas density by two orders of magnitudes does not 
change the inflow rate curves, however, the arm versus inter-arm 
overdensity 8p changes the inflow rate as illustrated by the sold 
and dashed curves. 

over three orders of magnitudes distribution. The value 
that we have calculated here for the dense gas streaming, 
at 40 pc distance from the SMBH in NGC 1097 corre- 
sponds to M ~ 0.033 MEddj where the Eddington accre- 
tion rate is onto a SMBH mass of 1.2 x 10 8 M (ba sed on 
the central stellar velocity dispersion adopted from |Lewis| 
& Eraclcous 2006J and the M,-a relation of Tremaine et 
al. 2002). 

6. CONCLUSIONS 

We have presented a detailed kinematics analysis of 
the dense interstellar gas in the circumnuclcar region of 
the nearby Seyfert 1 galaxy NGC 1097. We used ALMA 
Band 7 high-resolution observations of the HCN(4-3) 
emission line, which is postulated to trace tt-h 2 ~ 10 8 
cm -3 at 40 K kinetic temperature. While visual signa- 
tures of a rotation pattern dominate the observed veloc- 
ity field, we have successfully applied Fourier decompo- 
sition of the velocity field to unveil the kinematic sig- 
natures of two prominent spiral arms (in morphology). 
The number of nuclear spiral arms in NGC 1097 has been 
subject to a controversy and in the past years it has been 
clear that the presence of strong obscuration by dust at 
the very centre of this galaxy has complicated matters. 
We argue here that using velocity information from the 
HCN(4-3) clarifies matters in an unprecedented way. 

We have found a striking agreement between the kine- 
matics of the HCN(4-3) and an analytic spiral model 
that we previously build using ionized gas kinematic data 



at similar spatial resolution. The new ALMA data con- 
firm that the spiral arms have pitch angle of 52 ± 4°, 
down to < 40 pc from the SMBH in NGC 1097. 

We note that some studies have found that the HCN 
molecule could arise from radiation and vibrational ex- 



( T 
li 



citation (]Lepp & Dal garno 1996 Imanishi et al. 2004 



Meier & Turner 2012]). The relative contribution of the 



radiation pumping could then decrease the gas density 
with which the HCN is associated, and we note that an 
assumed density decrease by two orders of magnitudes 
will yield h/R < 0.01. For a thin disc, we have then mea- 
sured the mass inside the 100 pc radius to be 3.5J^'o x 10 8 
M Q , and inside the 40 pc radius, 8.0±|;§ x 10 6 M . The 
derived mass at 40 pc agrees with the < 3.0 x 10 7 M , de- 
rived from t he observed line flux fol lowing the procedure 
described in Gao & Solomon (2004). 

We have used a constant arm versus inter-arm overden- 
sity of 10% and kinematically derived a dense gas inflow 
of 0.3 M yr" 1 at 100 pc and 0.09 M Q yr" 1 at 40 pc ra- 
dius. In combination with our previously derived values 
from the ionized and CO gas, we calculate a molecular 
and ionized gas infall of 0.6 M© yr _1 at 100 pc and ~ 0.2 
M Q yr _1 at 40 pc distance from the central SMBH of 
NGC 1097. This inflow corresponds to M ~ 0.066M E dd 
onto a black hole in NGC 1097 with a mass of 1.2 x 10 8 
Mq. From the current data, it is not clear how much, 
if any at all, of the observed HCN(4-3) is radiationally 
excited by the active nucleus. In the presence of radi- 
ational excitation, the dense gas scale height that we 
present here will be lower limits. Notwithstanding, the 
gas inflow rates remain unchanged. 
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